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Structural distortions in cuprate materials offer a microscopic origin for anisotropies in electron 
transport in the basal plane. Using a real-space Hartree-Fock approach, we consider the ground 
states of the anisotropic Hubbard (t x 7^ t v ) and t-J (t x 7^ t y , J x 7^ J y ) models. Symmetrical but 
inhomogeneous ( "polaronic" ) charge structures in the isotropic models are altered even by rather 
small anisotropies to one-dimensional, stripe-like features. We find two distinct types of stripe, 
namely uniformly filled, antiphase domain walls and non-uniform, half-filled, in-phase ones. We 
characterize their properties, energies and dependence on the model parameters, including filling 
and anisotropy in t (and J). We discuss the connections among these results, other theoretical 
studies and experimental observation. 



I. INTRODUCTION 

The issue of inhomogeneous charge and spin order in 
cuprate systems, or "stripe phases" , has become one of 
the most controversial issues in the debate over high- 
temperature superconductivity. Experimentally, the ev- 
idence for static stripe phases in a variety of materi- 
als is incontrovertible: neutron diffraction studics_of the 
Nd-doped LSCO system La2-z-yNd y Sr x Cu04ErEI show 
over the superconducting range of a; a charge modula- 
tion which has been interpreted as arising from a system 
of charged, alternating (1,0) and (0,1) domain walls sep- 
arating antiferromagnetic (AF) regions where no holes 
are present. A similar phenomenon has been reported 
more recently in very low-doped LSCO without Nd, al- 
though in this case the stripes are diagonal and the sys- 
tem insulatingou The variation of this inhomogeneous 
charge structure with hole doping x has been character- 
ized from the incommensurate magnetic diffraction-peaks 
due to the antiphase nature of the domain walls .Bu Al- 
though dynamical stripes (below) have been proposed as 
a pairing mechanism,0 also beyond dispute is that static 
stripe formation is inimical to superconductivity, as ob- 
served both from the individual, possibly coexisting, or- 
der parameters at fixed Nd content yp and more broadly 
from the evolution of T c with yn 

Somewhat less unanimity is achieved concerning the 
existence and role of dynamical stripes, by which is meant 
charge fluctuations of the same incommensurate type. In 
contrast to easly, proposals based on the shape of the 
Fermi surface,Eni3 such fluctuations are now the favored 
interpretation of the incommensurate peaks—observed 
by ineLastic neutron scattering in both LSCOOli3'Q and 
YBCOEjo systems. ■— .Nuclear Quadrupole Resonance 
(NQR) measurementscjiill suggest that charge-ordering 
features are ubiquitous, in materials both with and with- 
out Nd, on the short timescales probed by this technique. 
Indirect evidence for dynamical stripes has also been in- 
ferred from the suppression of T c on impurity doping, by 
the argument that pinning effects lead to static stripe 



formation (see Ref. [18[ and references therein). 

On the theoretical side, indications for stripe phases 
predate their experimental observation by some 6 
years. Initial results from the real-space Hartree-Fopk 
approach applied to a three-band Hubbard modeO 
were followed by similar studies in a selection of 
systems,Ej~Ej and- by proposals based on frustrated 
phase separation.c3 A dramatic increase in this activity 
since the discovery of stripes in nickelates and cuprates 
has seen charge inhomogeneities investigated further in 
Hubbard- and t- J- based models by the Hartree-Fock 
(HF) approximationjHIH HF with additional correlations 
and static phonons,E3 the slave-boson technique,E3 nu- 
merical DensitvJilatrix Renormalization Group (DMRG) 
simulations,EirE3 the Dynamical Mean-Field. Theory 
(DMFT),H and Exact Diagonalization (ED) £2 The re- 
sults of all of these studies imply that one-dimensional 
(Id) charge ordering is a leading instability in this class of 
models throughout the physical parameter regime. Fur- 
ther analyses using the. nonlinear sigma model, LJ cou- 
pled Luttinger liquids, £3 and Cluster Perturbation The- 
ory (CPT)e3 adopt this striped state at the outset and 
shed additional light on its microscopic nature. 

At the microscopic level, the relationship between 
structure and superconductivity has also been in- 
vestigated since before the first reports, of charge- 
inhomogeneous phases. Buchner et aln noted the 
structural transformation and loss of superconductiv- 
ity on increasing Nd doping in a series of samples 
La2_ 2; _j,Ndj ) Sr 2 ;Cu04, and were further able to relate 
these results to the bond lengths and angle distortions 
in the low-temperature tetragonal (LTT) phase. Such 
analyses have been extended within the LSCO system to 
cover a variety of substituentsj^j and also towards the 
use of dynamical probes to investigate the relation of 
structure and spin correlations (see for example Ref. [39]) . 
Of particular interest in connection with our analysis is 
the LTT structural transition also caused by Eu doping, 
and the Cu NQR studiescJ of this phase which indicate 
the presence of three inequivalent Cu sites depending on 
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the location of the atom within or outside charged stripe 
structures. The microscopic relationship between lattice 
structure and the stripe instability has also been stud- 
ied recently by detailed experiments focusing .ori bond 
lengths, or micro-strain, in the CuO polyhedra.cll 

In this work we begin with the experimental motivation 
provided - and quantified - by Ref. g that real cuprate 
systems are often anisotropic, and that increase of this 
anisotropy supports static stripe formation. Theoreti- 
cally, this suggests that appropriate starting models for 
the underlying physics of the cuprate planes would be 
anisotropic Hubbard (t x ^ t y ) or t-J (t x ^ t y , J x =^ J y ) 
modelsJiH To investigate the relationship between this 
anisotropy and Id, static charge order, we will adopt the 
most straightforward method known to provide initial, 
qualitative insight, and conduct real-space HF studies of 
anisotropic Hubbard and t-J models. We aim to illus- 
trate the effects of anisotropics on the static charge and 
spin configurations, and also to elucidate the energetic 
contributions (kinetic vs. magnetic) which drive the for- 
mation of inhomogeneous structures. The unrestricted 
HF technique is by nature suitable only for discussing 
the static properties, as it neglects quantum fluctuations 
which are essential for dynamical properties and restora- 
tion of broken symmetries. Thus we will not address the 
contentious issues of dynamical stripes or the connection 
to superconductivity. However, we hope to provide a 
suitable foundation on which to base further consider- 
ations, and will discuss below some more sophisticated 
approaches to the anisotropic models which may shed 
light on these questions, as well as on the roles of quan- 
tum fluctuations and lattice degrees of freedom which are 
neglected at this level. 

The structure of the paper is as follows. In section 
II we review the experimental observations of structural 
anisotropics in cuprate materials. In Sec. Ill we give a 
brief account of the real-space Hartree-Fock method. Sec. 
IV presents results obtained for the Hubbard model with 
variable hopping anisotropy, and a comprehensive inves- 
tigation of their evolution with the other intrinsic and 
extrinsic system parameters. In Sec. V we consider the 
t-J model, provide analogous results where possible for 
anisotropics in both t and J, and compare the two mod- 
els. Sec. VI contains a discussion of the results in the 
context of experiment, some consideration of other mate- 
rials, models and theoretical approaches, and a summary 
of our analysis. 

II. STRUCTURAL DISTORTIONS IN CUPRATES 

Before the first observations of stripes in any oxides, 
doping of LSCO by rare-earth elements was found to 
cause a structural transition from the, low-temperature 
orthorhombic (LTO) phase to LTT.Hcj Biichner et 
al. investigated the suppression of superconductivity 
by Nd-doping in LSCOjj finding a phase diagram for 



La2- 2; -yNdj,Sr :E Cu04 with the LTO phase (Fig. 1(a)) for 
y < 0.16, LTT (Fig. 1(b)) with a superconducting ground 
state in the doping range 0.16 < y < 5x — 0.55, and LTT 
but non-superconducting for y > bx — 0.55. For large 
x the high-temperature tetragonal (HTT) phase remains 
stable down to low temperature, and for small x with 
y > 0.16, the structure is Pccn. By combining trans- 
port measurements with structural refinement, the line 
y ~ 5x — 0.55 was found to correspond to a critical buck- 
ling angle of the CuOe octahedra in the LTT structure, 
$ c ~ 4°. A further key point in support of the close con- 
nection between charge order and lattice structure is the 
observation!!! that the orientation of the stripes alternates 
between (1,0) and (0,1) on proceeding along the c-axis, 
as does the rotation axis of the LTT structural distortion 
(Fig. 1(b)). It is clear that the electronic structure of the 
Cu02 plane is altered strongly by small changes in the 
lattice structures 

The LTO and LTT phases both involve a distortion 
of the flat Cu02 planes by rotation of the CuOe struc- 
tural units. For LTO this rotation takes place around a 
(1,1) diagonal of the Cu02 system, such that all the in- 
plane O atoms of the HTT structure are displaced from 
the plane. For LTT this rotation is around a (0,1) or 
(1,0) axis, such that only x-axis (or y-axis) O atoms are 
displaced. Structural analysis indicates that the funda- 
mental Cu-0 separation is not altered significantly by 
the distortion, which is taken up rather by the Cu-O-Cu 
bond angle, 2$. Inspection of Fig. 1 shows that (de- 
spite the crystallographic terminology, which is based on 
a unit cell of side y/2 larger than the Cu-Cu distance), 
this planar buckling affects the Cu-O-Cu bonds equally 
in x- and y-directions in the LTO phase, but unequally 
in LTT. 

This result provides a microscopic origin for in-plane 
anisotropics in the LSCO system. Calculation from first 
principles of the electron transfer and superexchangc 
parameters, t and J, indicate that in cuprate materi- 
als these deviate rather. stronglycS from the elementary 
orbital-overlap theoryO For a (1,0) LTT anisotropy, one 
expects an angular dependence given approximately bycll 

t y = t x \ C0S(7T - 2$)|, Jy — J x COS 2 (7T - 2$). (1) 

For a maximal distortion angle! $ = 5°, the relative 
anisotropies \t x — t y \/t x ~ 1.5% and \ J X — J y \/J x ~ 3% 
may appear rather small, but we note that even these 
small values result in the absolute differences \t x — t y \ and 
I Jx ~ Jy \ exceeding the superconducting T c . The qualita- 
tive picture thus emerges of a possible microscopic origin 
for the formation of Id structures (stripes) as the true 
ground state of the weakly anisotropic electronic system, 
and it is this possibility which we wish to investigate. 
Competition of such a static stripe state with a super- 
conducting one, which may be the true ground state of 
the 2d system, would be fully consistent with experimen- 
tal observation. We note in this connection that the is*jp§ 
of bulk coexistence of these two statesa remains openBca 
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FIG. 1. Schematic representations of the conformation of 
the CUO2 plane in LTO (a) and LTT (b) lattice structures. 
Small, black circles represent Cu atoms, which remain in the 
plane in both structures. Large, black circles represent O 
atoms which remain in the plane, white circles O atoms which 
are driven above the plane, and grey circles O atoms which 
are driven below it by the tilting distortion. 



An important feature of the Nd-doped LSCO system 
which we have discussed is that the LTO-LTT transition 
and the formation of static stripes, or charge order, ap- 
pear to be decoupled. For y = 0.4 and in the optimal 
doping range of x, the structural transition occurs at a 
temperature Tut rff-SOK., while charge order appears at 
some T c h < 40K.BO Although these temperature scales 
become closer for some smaller values of y and x, they 
show opposing trends as functions of x, and it_appears 
justified to treat the two phenomena separately.ca We will 
thus consider only the energetics of the electronic subsys- 
tem, studied on a rigidly distorted lattice. This is clearly 
a simplifying feature in comparison with, for example, 
highly-doped nickelate systems, where stripe formation 
and structural distortion appear to occur together as a 
single, first-order transition, and any energetic consider- 
ations would require a comparison between the gain in 
the electronic subsystem with the cost to the lattice. 

In closing this section, we note that another manifestly 
anisotropic cuprate system is YBCO, where the 6-axis 
chains in the Cu(l)-0 layer cause a significant difference 
in the structural parameters of the CuC>2 plane. The 
absence of static stripes in YBCO suggests that x- and 
y-axis transfer integrals may in fact be rather similar (to 



within some critical ratio), or that the bilayer nature may 
be important, or may reflect the fact that the resultant 
stripe structure would be the same in every plane, pos- 
sibly suffering additional penalties in Coulomb energy 
compared to the case of LTT LSCO, where the rotation 
axis and stripe direction alternate along the c-axis. 



III. REAL-SPACE HARTREE-FOCK 
APPROXIMATION 

This technique has been employed by many authors to 
gain insight into the qualitative physics of the Hubbard 
model, and we provide only a brief overview of the steps 
involved in the procedure After the initial results of 
Zaanen and Gunnarsson,E3 who found stable, open and 
closed domain waUs-Ln-a. three-band Hubbard model, sub- 
sequent analysesEjt3'Efl of Hubbard models have verified 
this property in the one-band case, considered the evo- 
lution of stripe features with the ratio U / 1 and with fill- 
ing, and clarified both the filling of domain walls and the 
phase relationship between magnetic regions.EZI Further 
refinements of the method have involved the, inclusion of 
local correlations and static phonon terms. E3 

The one-band Hubbard model is written as 



H 



E 



t v (c. 



t , 



H.c. 



(2) 



where n.i a = cj^c^, 77 denotes x, y, and on an anisotropic 
lattice t x ^ t y . Denoting the terms in Eq. (||) by H t and 
Hjj , in the HF approximation 



H v = U^2 [{ n ij) n H + "-it( n a> ~ ( n il)( n n) ( 3 ) 

i 



where the signs include fermion statistics, precluding the 
need to form Slater determinants of quasiparticle states. 
We will focus here on collinear .scin configurations, and 
not in spiral or vortex solutions.cil In this case, by choos- 
ing the axisj-of the possibly broken spin symmetry to be 
the z- axis,E3 one may work without loss of generality 
using only the first line of Eq. (||) . 

The model is solved in real space on a small cluster 
of sites i. The particle number qi and magnetization 
components mf are defined by the expectation values 



Qi = 


( c !t c *t + c !i c a) 


(4a) 


mf = 


2 ( c !| c jT — 4l c il) 


(4b) 


m x i = 




(4c) 


m y L — 


-^( c It c 'U - c ^ T l Ci T>• 


(4d) 



Solution proceeds by iteration of the site parameters to 
self-consistency, at which point one has solved all of the 
mean- field equations contained in Eq. (0), subject to the 
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FIG. 2. Numbering scheme for an / = n x n cluster, 
illustrated with n = 4. 



additional constraint on the total particle number, which 
is implemented in the number of filled quasiparticle en- 
ergy levels. This procedure is not identical to finding 
the global minimum of the multi-parameter free-energy 
surface, in that in principle it may only find stationary 
points, and the user is generally offered no guarantee that 
these correspond to global minima. 

To solve the self-consistency problem for an I = n X n 
cluster requires diagonalization of the 2n 2 x 2n 2 Hamil- 
tonian matrix form of Eq. (0) , 



Ni Mi 2 M*! Mis 



M? 2 



N 2 
M 
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M 23 
N 3 M 34 



(5) 



Here we have discarded the constant terms in Eq. (|3|) , al- 
though these are restored in calculating the energy of the 
self-consistent solutions which emerge. The 2x2 "diago- 
nal" and "off-diagonal" matrix elements are respectively 



N; 



\U{q t - 2m?) 
—U(rrif + imf) 



— U(m\ 



- im; 
2m? 



and 



iVliJ "lot, 



(0) 



(7) 



in which rj = x,y for j = i + 1(— n), i + n(—l), and the 
two-component structure corresponds to the two spin di- 
rections in the S z basis. For the numbering scheme il- 
lustrated in Fig. 2 with a 4 x 4 cluster, the off-diagonal 
blocks connect i to i ± 1 in the ^-direction, and to i ± n 
in the y-direction. Open boundary conditions (BCs) are 
implemented with no forward (backward) connections in 
x when mod(i, n) = (mod(i, n) = n— 1), and in y when 
i > n 2 — n (i < n). Periodic BCs have forward (back- 
ward) connections from i to i + 1 — n (i — 1 + n) in i, 
such as the element M41 in Eq. (||), and to i + n — n 2 
(i — n + n 2 ) in y for the same edge sites. Because the 
Hamiltonian matrix is Hermitian, only half of these con- 
nections must be set explicitly. For cluster filling x, the 
lowest 2n 2 x eigenstates of the diagonal Hamiltonian ma- 
trix are filled, g, and for all sites i are deduced from 
the eigenvectors representing the new quasiparticles, the 
new matrix is diagonalized, and the process continued. 
The consistency condition was generally taken to have 



been achieved when the summed squares of the changes 
in <7j and mf were less than 10~ 8 per site, although this 
was continued in some cases to the machine precision 
to ensure that no further changes were possible. Conver- 
gence was usually achieved for 12 x 12 systems within 200 
iterations using open BCs, and within 1000 iterations for 
periodic BCs, where the iteration procedure could more 
clearly be observed converging on and diverging from a 
sequence of local minima. 

The t-J model has previously been considered in the 
real-space HF formalism only in Ref. |2^. Here we have 
chosen to study the form 



H t -j = H t + Hu + J (Sj-Sj- - jUiUj) , 



(8) 



with rii = ritf + riii, anc ^ ^ taken to have a large, fi- 
nite value to mimic the projection property of the t-J 
model onto only those states with no doubly occupied 
sites. We will have more to say on this approximation 
in Sec. V. In this case the J term admits further de- 
couplings, among which we neglect pairing terms when 
working in a canonical ensemble, and nearest-neighbor 
spin-flip terms because these are generally small. Thus 
we retain only a finite expectation value of the bond hop- 
ping term A ija = (c\ a Cj a )m 

The HF decomposition of the additional J term is 



Hi = -U 



1 tVM 



+^-TsVt + A h c u c 3i + A h A ai ( 9 ) 

+ A % j t c ] x c * T + A ij t c li c o T + A* 3 T Aij 1 . 
The 2n 2 x 2n 2 matrix to be diagonalized (cf. (||)) has 

/ ±{U( qi -2mf)- J lL ] -U(mf-imy) \ 
1 \ -U{mx + im v t ) l[U( qi + 2ml)- J lT ] J 

(10) 

as the 2x2 diagonal component, where 

J o Jx(j^i-\-xa ^i — xa) "T" Jyij^i+ycr T^i— yo)i (^^-) 

and 





13 \ n t - i j A* 

u L n 2"''/ i)T 



(12) 



in the non-zero, off-diagonal blocks of the upper triangle. 
Further details of the iterative calculations in this case 
are deferred to Sec. V. 



IV. HUBBARD MODEL 

In this section we describe the charge and spin config- 
urations (|) which emerge as the self-consistent solutions 
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to the Hubbard model (g) in the Hartree-Fock approxi- 
mation (H). We begin by considering the isotropic system 
as a function of the "intrinsic" system parameters, the ra- 
tio U/t, the (hole) filling x, and the temperature T, and 
make contact with the results of previous studies. We 
then illustrate how these results are altered in the pres- 
ence of hopping anisotropy t x ^ t y , for which we study 
both the charge and spin distribution and the kinetic and 
magnetic energy components. We comment also on the 
variation of the solutions with "extrinsic" parameters, by 
which is meant cluster size, BCs, and commensuration. 
We now set the energy unit to be t = 1. In all figures to 
follow the i-axis is horizontal and the y-axis vertical. 



A. Variation of U 

For small values of U the solutions are homogeneous in 
charge and have no local magnetization. This "metallic" 
phase is the ground state for [7 < 3, the exact value of 
U c depending on the filling. We have not analyzed this 
regime or the metal-insulator transitioH-.in any detail, 
as many more accurate techniques exist .Ej For all values 
U > 3, the predominant HF solutions are inhomogeneous 
charge clusters, with locally AF spin orientation. In the 
approximate range 3 < U < 6, these clusters take the 
form of closed loops of electron-depleted sites, or domain 
walls enclosing an undoped region (Figs. 3(a) and 4(a)). 
We will call these formations "corrals", to distinguish 
them from the situation as U is increased. For U > 6 the 
charge clusters are better described as a solid region of 
hole doping, which we will call a "polaron" (Fig. 3(b)). 
Our results are rather similar to those in Fig. 1 of Verges 
et a/.Ea for this isotropic situation. 



B. Variation of x 

As the filling is varied in the isotropic model we 
note the same general features, but with important 
alterations. With smaller hole doping the cluster is 
less influenced by the system boundaries, and prefers 
(Fig. 4(a)) "diagonal" corrals composed of (1,1) domain 
walls, rather than the "vertical" type with (1,0) do- 
main walls (Fig. 3(a)). This result is quite ubiquitous, 
and readily understood from the kinetic energy gain due 
to nearest-neighbor hopping processes on all four bonds 
around one hole. Furthermore, applying periodic rather 
than open BCs leads to a division of large corrals into 
small, diagonal ones, which appear to form a lattice 
(Figs. 4(b), 5(a)). We may conclude that the vertical 
corral shown in Fig. 3(a) is constrained by the system di- 
mensions, an observation confirmed at still higher filling. 
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FIG. 3. Ground-state charge distributions for Hubbard 
model with hole doping a; = l/6onal2xl2 cluster with open 
BCs. (a) U = 5, t x = t v = -1. (b) U = 8, t x = t y = -1. 
(c) U = 5, t x = -1.05, t v = -0.95. (d) U = 8, t x = -1.05, 

-1.1, 



(e) U = 5, t x 



t U 



-0.9. (f) U = 8, 



tx — 



-1.1, t y = —0.9. In these figures the largest circles 



correspond to (rij) = 0.48, or 52% hole doping of the site. 
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FIG. 4. Ground-state charge distribution for Hubbard 
model with hole doping x = 1/8 and 17 = 5 on a 12 x 12 
cluster, (a) Open BCs, = t y — — 1. (b) Periodic BCs, 
t x = t y = -1. (c) Open BCs, = -1.05, t„ = -0.95. 
(d) Periodic BCs, t a = -1.05, t» = -0.95. (e) Open BCs, 



t*E 1-1, tu 



Charge scale as in Fig. 3. 



-0.9. (f) Periodic BCs, t x 



-1.1, t v 
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C. Variation of temperature 

We have in addition considered the real-space HF ap- 
proximation at finite temperature. This part of the anal- 
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FIG. 5. 
models on 
odic BCs, 
U = -1.1 



Ground-state charge distribution for Hubbard 
a 12 x 12 cluster, (a) U = 5, x = 1/6, peri- 

t x = t y = -1. (b) U = 8, x = 1/8, periodic BCs, 

t y = —0.9. Charge scale as in Fig. 3. 



ysis was motivated by the picture of the self-consistency 
procedure as the following of a path through a free-energy 
"landscape". Because the solution is a stationary point, 
and not a global minimum, one would wish to sample as 
many local minima as possible, which is facilitated when 
the landscape be as smooth as possible. On descending 
from high temperatures, one expects a uniform charge 
distribution to be the state of minimum free energy, and 
that the leading instability to an inhomogeneous distribu- 
tion would appear as a single minimum before the land- 
scape could become more complex at lower temperatures. 
The effect of a finite temperature is readily incorporated 
by a thermal factor in each of the mean-field equations 
(Q), and by replacing the condition on the number of 
quasiparticle states filled with an equation determining 
the Fermi energy. 

At high temperatures, T > 2, the solutions are al- 
ways uniform in charge and have no magnetization. De- 
pending on the exact choice of the other parameters, for 
temperatures in the range 0.18 < T < 2, an AF spin con- 
figuration develops while the charge distribution remains 
uniform. Finally, at temperatures T ~ 0.18 for most 
of the situations we have considered, an inhomogeneous 
charge structure develops (corrals, polarons, stripes with 
suitable anisotropy) while the spins remain locally AF. 
In all the cases we considered, the leading instability to 
charge inhomogeneity at finite temperature was identical 
to the ground-state structure found by working entirely 
at T = 0, with a single exception mentioned below. This 
result is reassuring in the sense of helping to confirm that 
the minima emerging from the zero-temperature studies 
are quite likely to be global in nature, and in this work 
we will not pursue the finite-temperature studies beyond 
establishing this point. We note finally that in this ap- 
proach the onset temperature for formation of an incom- 
mensurate spin structure where this arises (e.g. Fig. 3(e)) 



is identical to that for charge order. 

D. Variation of anisotropy 

We turn now to the most important aspect of our anal- 
ysis. Figs. 3 and 4 (c)-(f) illustrate the effects of in- 
creasing the hopping anisotropy e t = t x /t y — 1 to 11% 
and 22% for a selection of initial parameters and BCs. 
Figs. 3(a), (c) and (e) give the clearest demonstration 
of stripe formation as a result of hopping anisotropy; al- 
though this is for open BCs the results for periodic BCs 
are almost identical for this choice of U and x. Figs. 3(b), 
(d) and (f) show that for stronger U the polaronic state 
is more favorable, and although it becomes progressively 
more elongated with increasing anisotropy, stripe forma- 
tion is not achieved at 22%. Similarly, Figs. 4(a), (c) 
and (e) show the tendency of corrals to become increas- 
ingly diamond-shaped with anisotropy, before crossing 
to a fully Id state. The difference between these and the 
analogous cases in Fig. 3 is that for the smaller filling the 
(diagonal) corrals are rather more stable, and a greater 
anisotropy would be required to create the stripe state. 

One of our key qualitative results is evident in these 
figures. The direction of the stripes, and the major axis 
of the diamonds, is perpendicular to the direction of the 
strong hopping. Stripes are stabilized not by charge mo- 
tion along them (the pictu*G.-af a conducting channel), 
but by transverse hopping,E3l23 onto and off the stripe. 
This is readily shown by a simple argument based on site 
charge densities, which we illustrate for the 4-site clus- 
ter with two different densities n\ and n 2 of Fig. 6. For 
both spin species (nj ,n\ = n), a hopping amplitude pro- 
portional to the occupation of the initial state and the 
availability of the empty state gives 

(££)/* = ni(l-ni)+na(l-n2) 

= 2(m + n 2 - (nj + nj)) (13a) 
(E t v )/t = ni(l-n 2 )+n 2 (l-n 1 ) 

= 2(ni+n 2 -2mna), (13b) 

whence the identity nl + n^ > 2nin 2 ensures that (E*) > 
(El.). That hopping between two different charge densi- 
ties is favored over that between equal densities ensures 
that transverse charge fluctuations determine the orien- 
tation of uniform stripe states such as those in Fig. 3(e). 

This result should not be confused with the metallic 
or insulating nature of the stripes: these are still sep- 
arated in the transverse direction by insulating, AF re- 
gions. Of all the charge-inhomogeneous solutions illus- 
trated here, only the stripe solutions are metallic, and 
this in the stripe direction within the 2d structure. While 
we have considered transverse hopping as uncorrelated, 
virtual processes, other microscopic approaches suggest 
a coherent hopping of stripe segments of lengths around 
5 lattice sitesH 



y i O O n, 

FIG. 6. Four-site cluster with two separate charge densi- 
ties ni and 712, to illustrate kinetic energy difference between 
x- and y-directions (Eq. (p^)). 

The above reasoning clarifies the somewhat unex- 
pected result in Figs. 4(d), 4(f) and 5(b): here the stripes 
are aligned in the direction of the stronger hopping. In 
this case the hole count at x = 1/8 is insufficient for uni- 
form stripe formation on a 12 x 12 cluster, and instead 
we find a non-uniform charge configuration. Such alter- 
nating stripes were presented in Fig. 1(c) of Ref. ^9] (to 
be contrasted with the variant in Fig. 1(a) of the same 
work), and are one HF analogue of "half-filled" stripes.cll 
We will follow this terminology, but note in passing that 
the idealized stripe with half of its sites occupied by holes 
is in fact quarter-filled with electrons. The fact that hop- 
ping along the stripe proceeds between sites of high and 
low charge density helps to explain ([l3]) the orientation 
of this type of Id charge inhomogeneity. 

An enduring problem for HF studies of the Hubbard 
model has been that these show a definite preference for 
uniform, filled stripes, as in Fig. 3(e) where the hole den- 
sity per site along the stripe is unity. By contrast, the ex- 
perimental situation favgrs half-filled stripes, in the inter- 
mediate doping range. Era While there is to date no infor- 
mation to distinguish between the possibilities of uniform 
or alternating (Fig. 4(b)) hole distribution within such a 
stripe, one may expect th&Jatter to be favoxed by increas- 
ing U. Some analytica£j£J and numericalEil studies have 
addressed the issue of the additional physics which may 
be responsible for this result. Here we have found that 
hopping anisotropy represents an additional factor which 
may contribute to the stability of non-uniform, half-filled 
stripes. Although this phase is obtained at the "magic" 
filling x = 1/8, close to which charge-ordering effects in 
experiment are strongest JI3 we hesitate to comment on 
a connection due to the issue of commensuration effects 
on the finite cluster. 

In common with other HF studies, this state is some- 
what more fragile than the uniform stripe, as indicated by 
the fact that it is preferred over a corral only in the pres- 
ence of periodic BCs, and has a smaller energy gain per 
hole than the uniform stripe (Table III below) . This case 
(Fig. 4(f)) was the only situation in which the leading 
instability at high temperatures (a corral) was different 
from the low-T structure. However, our result remains 
indicative of the possibilities allowed in the presence of 
hopping anisotropies, and we believe representative of the 
situation for appropriate fillings. At larger values of U 
one observes the emergence of the same type of state, but 
also (Fig. 5(b)) that competition is more severe from a 
phase of isolated, small polarons (Fig. 1(d), Ref. ^9j). In- 
creasing anisotropy promotes alignment of these polarons 
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(b) 

FIG. 7. Ground-state spin configuration for Hubbard mod- 
els on a 12 x 12 cluster with [7 = 5 and open BCs. (a) x = 1/8, 

= tj/ = — 1; corresponding charge distribution in Fig. 4(a). 
(b) x = 1/6, tc = —0.9, t y = —1.1; corresponding charge dis- 
tribution in Fig. 3(e). Length of largest spins represents 85% 
of full moment. 



towards the non-uniform stripe phase. 

A further qualitative point concerns the magnetic na- 
ture of the domain walls. As illustrated in Fig. 7, all of 
the corral and uniform stripe phases are found to have 
antiphase domain walls. Those sites located in the walls 
possess no magnetic moment, and the AEpXegimes on 
each side have a mutual phase shift of 7r£3'E3 Fig. 7(a) 
shows the spin configuration for the diagonal corral state 
in Fig. 4(a), illustrating from the phase shift between the 
inside and outside of the corral that the diagonal domain 
wall, whose sites have no moment, can be considered to 
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FIG. 8. Ground-state spin configuration for Hubbard 
model on a 12 x 12 cluster with U = 5, a; = 1/8, periodic 
BCs and t x = —0.9, t y — —1.1; corresponding charge distri- 
bution in Fig. 4(a). Spin scale as in Fig. 7. 



be filled, or uniform. Fig. 7(b) shows that the uniform 
stripe state in Fig. 3(e) possesses the same property. By 
contrast, we see in Fig. 8 for the non- uniform stripe state 
close to half- filling that there is no phase change (cf. 
Fig. 1(c), Re{Up9|), meaning that the stripe is an in-phase 
domain wall.E°TThis result is also opposite to that of the 
unifpjjmly half-filled stripe, which is an antiphase domain 
wallcJ'Eil for the same reason as for the filled stripe. In 
the non-uniform case, the sites in the center of the stripe 
retain a magnetic moment, and these are ferromagneti- 
cally (FM) aligned. However, we believe this to be driven 
not by the preference for a FM hopping channel, but by 
the spin orientations of the sites surrounding those with 
higher hole densities: these sites are analogous to a miss- 
ing electron, with charge and spin which reflect the hop- 
ping from the four nearest neighbors. The one-particle 
physics suggested by these observations, supports an in- 
terpretation by Nagaoka's theorem.EjO It appears that 
the non-uniform stripe could be considered as aligned, 
small FM polarons. 

Within the limits of the HF approximation, we may 
quantify some of the above observations by inspecting 
the energies of the solutions. Table I illustrates the 
effect of increasing anisotropy on the total energy for 
one parameter set (Figs. 3(a), (c) and (e)) by decom- 
position into potential and kinetic energy components. 
We see that increasing anisotropy causes a small in- 
crease in the potential contribution (double-occupancy 
energy), but that this is more than compensated by gains 
in the kinetic part. The strong effect of the hopping 
anisotropy is clearly demonstrated by the fact that for a 
22% anisotropy in t values (right column), the kinetic en- 
ergies differ by almost a factor of 2. This strong increase 
in energy difference reflects the dominance of transverse 
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TABLE I. Total ground-state energy E to t, and kinetic 
(Ex, E V K ) and potential {Ep) energy components per site for 
Hubbard model on a 12 x 12 cluster with U = 5, x = 1/6, open 
BCs and increasing anisotropy t x /t y ; corresponding charge 
distributions in Fig. 3(a), (c) and (e). 





tx j ty 1 


1.05/0.95 


1.1/0.9 


E^t 


-0.78782 


-0.79361 


-0.80317 


Er 


-1.22931 


-1.24634 


-1.25942 


E X K 


-0.61465 


-0.72917 


-0.80513 


El 


-0.61465 


-0.51717 


-0.45429 


Ep 


0.44148 


0.45273 


0.45625 



hopping in stabilizing the uniform stripe phase. 

Table II shows the total energies per site for 12 x 12 
clusters at hole dopings x of 0, 1/8 and 1/6, and Table III 
the differences between doped and undoped systems pre- 
sented as energy per hole. From Table II, systems with 
periodic BCs have significantly lower total energy, but a 
rather similar energy per hole. The same evolution of 
the total energy with anisotropy as in Table I occurs for 
all parameter sets with open BCs, and also at half-filling 
(x = 0), for the same reason as above. However, with 
periodic BCs a small anisotropy appears in most cases 
to result in a higher overall energy, but larger anisotropy 
values restore the behavior observed with open BCs. This 
result we ascribe to the competition of different possible 
structures within the system before a truly Id state may 
be established; periodic BCs vastly increase the degener- 
acy of available solutions, sometimes leading to conver- 
gence difficulties, and one has less confidence of finding 
the true ground state. Turning to Table III, the value of 
the energy per hole reflects both the mobility of the holes 
and the gain in potential energy due to a lower proba- 
bility of double occupancy. These trends compete as U 
is increased, and also as x is increased, as a consequence 
of which the values in the table vary rather little; this 
statement applies also for smaller (x = 1/12) and larger 
(x = 1/4) doping levels (not shown). 

In concluding this subsection, we have found that hop- 
ping anisotropy in a Hubbard model can give rise to two 
qualitatively different types of stripe state, whose appear- 
ance depends on the intrinsic system parameters, and 
whose nature is influenced by the commensurability of 
the filling with the cluster size. Uniform stripes have 
the same charge density and no magnetic moment on ev- 
ery site, align perpendicular to the direction of strong 
hopping, represent antiphase domain walls between AF 
regions on either side, and are energetically rather sta- 
ble. Non-uniform stripes have alternating charge density 
and finite magnetic moment on the central sites, align 
in the direction of strong hopping, and act as in-phasc 
domain walls. Both types of stripe are stabilized by hop- 
ping between sites of differing charge density, which in 
the uniform case means transverse hole hopping, but in 
the non-uniform case results in similar contributions from 



TABLE II. Ground-state energies per site for Hubbard 
model on a 12 x 12 cluster. O and P denote open and pe- 
riodic BCs, anisotropics quoted as ratio t x /t y . 



BCs 


u 


tx /ty 


x = 


1/8 


1/6 


O 


5 


1.00 


-0.6321 


-0.7514 


-0.7878 


O 


5 


1.11 


-0.6340 


-0.7534 


-0.7936 


O 


5 


1.22 


-0.6397 


-0.7563 


-0.8032 


P 


5 


1.00 


-0.6820 


-0.7967 


-0.8325 


P 


5 


1.11 


-0.6840 


-0.7905 


-0.8330 


P 


5 


1.22 


-0.6902 


-0.7950 


-0.8551 


O 


8 


1.00 


-0.4294 


-0.5550 


-0.5999 


O 


8 


1.11 


-0.4306 


-0.5591 


-0.6008 


o 


8 


1.22 


-0.4341 


-0.5624 


-0.6113 


p 


8 


1.00 


-0.4659 


-0.5974 


-0.6429 


p 


8 


1.11 


-0.4672 


-0.5971 


-0.6401 


p 


8 


1.22 


-0.4710 


-0.6002 


-0.6453 



TABLE III. Ground-state energies per doped hole for Hub- 
bard model on a 12 x 12 cluster at fillings x = 1/8 and x = 1/6. 
O and P denote open and periodic BCs, anisotropics quoted 
as ratio t x /t y . 



BCs 


u 


t X /ty 


£ = 1/8 


x = 1/6 


O 


5 


1.00 


-0.9539 


-0.9342 


O 


5 


1.11 


-0.9553 


-0.9575 


O 


5 


1.22 


-0.9328 


-0.9809 


P 


5 


1.00 


-0.9176 


-0.9032 


P 


5 


1.11 


-0.8502 


-0.8938 


P 


5 


1.22 


-0.8385 


-0.9890 


O 


8 


1.00 


-1.0052 


-1.0228 


O 


8 


1.11 


-1.0284 


-0.9909 


O 


8 


1.22 


-1.0268 


-1.0634 


P 


8 


1.00 


-1.0523 


-1.0619 


P 


8 


1.11 


-1.0391 


-1.0376 


P 


8 


1.22 


-1.0336 


-1.0461 



both longitudinal and transverse hole motion. 



E. Variation of extrinsic parameters 

In the preceeding subsections we have already noted 
some of the effects of the extrinsic parameters, a term we 
use to refer to cluster dimensions, BCs, and commensu- 
ration between filling and cluster size. Open BCs cause 
a definite tendency for holes to avoid the edges, where 
their hopping energy is lowered. While for small dop- 
ing this merely sets the inhomogeneous charge structure 
in the center of the cluster, at larger dopings it can act 
to influence the alignment of domain walls (Fig. 3(a)), 
particularly at smaller values of U where the charge 
structures are extended. Periodic BCs relax this con- 
straint, and thus facilitate the formation of domain lines 
with the optimal orientation (diagonal (Fig. 5(a)) or 
vertical (Fig. 4(d), (f), also as in Fig. 3(e))) for the 
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(a) (b) 

FIG. 9. Ground-state charge distribution for Hubbard 
models on a 12 x 6 cluster with U = 5 and t x = t y = — 1. (a) 
a; = 1/6, open BCs. (b) x = 1/8, periodic BCs. Charge scale 
as in Fig. 3. 

anisotropy, and clusters of the optimal size and configu- 
ration (Figs. 4(b), 5(a)). These effects are also insepa- 
rable from the issue of the commensuration between the 
cluster size and the number of doped holes, which is in- 
strumental in determining the type of the polaron lat- 
tice (Figs. 4(b), 5(a)), and the filling and nature of the 
stripes (Figs. 3(e), 4(f)), in the cases where these form. 
As in all such studies on finite clusters, the possibility 
arises of large, finite-size gaps in the quasiparticle spec- 
trum (for example, between fillings of 24 and 25 holes 
in the 12 x 12 system). However, we do not expect that 
such effects would particularly favor one class of charge- 
inhomogeneous solution over another. 

A more primitive effect of cluster size emerges when at 
least one of the dimensions becomes small on the scale of 
the equilibrium, inhomogeneous charge structure. This 
is illustrated in Fig. 9(a) for a 12 x 6 system with hole 
doping x = 1/6. Although the intrinsic parameters are 
the same as in Fig. 3(a), the inhomogeneous charge struc- 
ture of minimum energy in this system is a pair of an- 
tiphase domain walls, even for open BCs. This obser- 
vation should be borne in mind when interpreting the 
results of small-system studies such as exact diagonal- 
ization and DMRG, in which rather more sophisticated 
analyses are necessary to establish the true ground state 
of a hypothetical infinite system. Fig. 9(b) shows the 
analogous result for x = 1/8, where the 9 holes form one 
and a half stripes with periodic BCs. 

Finally, the nature of the ground state is very strongly 
influenced by the choice of periodic or antiperiodic BCs. 
This is more readily implemented by retaining periodic 
BCs, but working on a cluster with one side chosen 
to be odd, .and was the method by which Zaanen and 
GunnarssonE3 first observed stripe solutions in the form 
of a single domain wall. The analogous result for sys- 
tems even on both sides was a closed domain wall like 



the corrals presented above, but the straight domain line 
was found to offer the minimum energy per particle. We 
have also considered even x odd and odd x odd systems 
with both periodic and mixed BCs (open in the odd di- 
rection, periodic in the evened), and in the same way find 
the possibility of "trapping" stripe solutions in the form 
of single domain walls. Such studies are useful in the 
context of investigating the effects of anisotropy, as one 
may monitor the total energy of the fixed configuration 
as a function of t x /t y , and we will employ this approach 
below for the t-J model. However, for a variety of rea- 
sons we do not consider single, infinite domain walls to 
be viable solutions for the isotropic system in the ther- 
modynamic limit, and so have chosen not to dwell upon 
this type of analysis. 

For any 3d system of weakly coupled planes, long- 
rangerAF order is favored if the domain walls form closed 
loops S3 We note in our studies a preference for domain 
lines in corrals to lie diagonally, in accordance with ex- 
pectation for nearest-neighbor hopping, but that in no 
cases did a diagonal stripe phase appear with periodic 
BCs. Instead we find a preference for small, diagonal po- 
larons rather than extended lines in the isotropic system, 
and that vertical stripes required the LTT-inspired hop- 
ping anisotropy. A further qualitative observation in this 
direction is the absence in our isotropic, evenxeven stud- 
ies of a phase with two parallel, vertical domain walls, 
even at the optimal filling for this state (cf. Fig. 3(e)). 
With the goal of elucidating the effects of anisotropy on 
the static solutions, we have concentrated on the most 
unbiased initial conditions which appear to retain the 
possibility of sampling the candidate ground states. 

V. t-J MODEL 

We turn now to anisotropy effects in the t-J model (j^). 
The previous section has given considerable insight into 
the complex range of possibilities which may arise. Nu- 
merical studies of the t-J mode£2l have provided evidence 
for charge structures of closed domain-wall loops similar 
to those we find in the Hubbard model, but also for a 
close competition with many other structures, including 
open lines (stripes). In order to assess anisotropy effects, 
and to address the question of collaboration or competi- 
tion between t and J in stripe formation and alignment, it 
will be necessary to characterize the favorability in differ- 
ent regimes of 2d vs. Id structures, corrals vs. polarons, 
in-phase vs. antiphase domain lines, and the possible role 
of phase separation.E3~E£l 

In the real-space HF approximation to the t-J model 
one encounters an immediate technical problem. The 
value of U shold be very large (Sec. Ill), in order to 
have the Hubbard term act to project out any states 
with doubly-occupied sites from the space of those con- 
sidered. In a selection of studies with U = 100, we found 
the results to vary appreciably with changes in J, and 
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(a) (b) 

FIG. 10. Ground-state charge distribution for Hubbard 
models on a 12 x 12 cluster with x = 1/6, open BCs, and 
t x =t y = -l. (a) U = 100. (b) U = 1000. Charge scale as in 
Fig. 3. Spin distributions (not shown) are AF with unaltered 
phase. 

particularly with changes of J x relative to J y , on the or- 
der of 0.05. The "residual" superexchange term in a pure 
Hubbard model, J = 4t 2 /U = 0.04 for this choice of U, 
suggests that still larger values of U are required for a 
degree of projection which can safely be taken to repre- 
sent the asymptotic limit. This observation is supported 
by computing the value of the double-occupation energy. 
Indeed, as shown in Fig. 10, while a Hubbard model with 
U = 100 gives a robust, polaronic state with variations 
of charge between sites, the situation for U = 1000 (and 
larger) corresponds to a cluster of uniformly charged sites 
more reminiscent of incipient phase separation, and with 
no obvious contribution from kinetic terms. We note 
also in connection with Fig. 10(b) that the calculation 
has converged to a state lacking the fourfold symmetry 
expected from a uniform initial configuration, which is 
indicative of difficulties in employing the real-space HF 
technique in this parameter regime. 

For such large values of U, (i) every local minimum 
becomes very deep (0(U)), and (ii) the fluctuation ener- 
gies (0(t 2 /U)) become very small. Calculations in this 
regime show that the method loses its flexibility to find 
a global minimum, and that the final state is determined 
largely by the chosen starting configuration. Physically, 
large U projects out the relevance of quantum fluctu- 
ations, a statement we will make more specific in the 
following example. There remain two possible courses of 
action by which one may gain insight into the anisotropic 
t-J model using the real-space HF approach. As outlined 
in the preceeding section, the first is to create a captive 
stripe at large U, and to study through the ground-state 
energy E to t the qualitative question of how anisotropics 
in t and J affect such a state. The second is to reduce 
the value of U to a size at which the system appears once 
again to have the ability to explore a number of distinct 



(a) (b) 

FIG. 11. Ground-state charge distribution for t-J models 
with periodic BCs, t x = t v = —1, J x = J v = 0.35, and 
U = 1000. (a) 12 x 11 cluster with 12 holes, circle radius on 
stripe represents 98.1% hole occupation; (b) 12 x 12 cluster 
with 6 holes, circle radii on stripe represent 55.9% and 20.3% 
hole occupation. 

states, and to seek indications for the validity of the re- 
sults in a model with true projection. 

A. Large U 

In this subsection we consider the kinetic and magnetic 
energies of filled and half-filled stripes as a function of the 
anisotropy. We set U — 1000 to ensure adequate removal 
of doubly occupied sites, and consider the canonical t- 
J model parameter ratio J/t = 0.35. We trap stripes 
of 12 and 6 holes along the center line of a 12 x 11 and 
12 x 12 clusters, respectively, with periodic BCs, and com- 
pute the kinetic and magnetic energies of the system with 
±11% anisotropy in both t and J. Because the stripe is 
aligned along the ^-direction, automatically breaking any 
4-fold symmetry, positive and negative anisotropy effects 
are not necessarily the same. We find that the schematic 
picture of the fully local stripe is very nearly realised 
(to 98%, Fig. 11(a)) in the uniform case (12 holes) when 
the spin configuration is set such that the stripe is an 
antiphase domain wall, as in the previous section. Sim- 
ilarly, the half-filled stripe solution (Fig. 11(b)) emerges 
only when the stripe is an in-phase domain wall, which 
accounts for the different transverse cluster dimensions 
used in each case. Thus we do not show the correspond- 
ing spin distributions for Fig. 11, which consist simply of 
full-moment AF regions on both sides of the stripe, with 
(a) and without (b) a phase shift across the domain wall. 

The results for kinetic and magnetic energies are pre- 
sented in Tables IV and V. We list these separately be- 
cause the numbers are not directly comparable at large 
U : while the magnetic energy is determined essentially 
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TABLE IV. Kinetic (upper figure) and magnetic (lower) 
energy components per site for t-J model on a 12 x 11 cluster 
with 12 holes, periodic BCs, and U = 1000. 





t x — —0.95 
t y = -1.05 


t x = -1 
t y = -1 


t x = —1.05 
t y = -0.95 


J x = 0.3325 
J y = 0.3675 


-0.045534 
-0.30086 


-0.043226 
-0.30088 


-0.040996 
-0.30090 


Jx = 0.35 
Jy = 0.35 


-0.045512 
-0.30163 


-0.043206 
-0.30165 


-0.040976 
-0.30167 


J x = 0.3675 
J y = 0.3325 


-0.045491 
-0.30239 


-0.043185 
-0.30241 


-0.040955 
-0.30243 



TABLE V. Kinetic (upper figure) and magnetic (lower) en- 
ergy components per site for t-J model on a 12 x 11 cluster 
with 6 holes, periodic BCs, and U = 1000. 





tx = -0.95 
t y = -1.05 


tx = -1 
ty = -1 


t x = -1-05 
t y = -0.95 


J x = 0.3325 
Jy = 0.3675 


-0.089029 
-0.31249 


-0.088635 
-0.31261 


-0.088445 
-0.31272 


Jx = 0.35 
Jy = 0.35 


-0.089018 
-0.31255 


-0.088626 
-0.31267 


-0.088438 
-0.31280 


Jx = 0.3675 
J y = 0.3325 


-0.089006 
-0.31260 


-0.088615 
-0.31275 


-0.088429 
-0.31289 



by J, and depends only on the number of antiparal- 
lel spin pairs, the kinetic energy requires a transfer of 
electrons between sites. Although a transverse hopping 
would appear to be favorable between neighboring filled 
and empty sites, the expectation value {c irT Cj a ) depends 
on the component of the final-state configuration present 
in the HF wave function (which should be considered as a 
linear superposition of many charge configurations in the 
basis of lattice sites) . Large values of U act to suppress 
these additional components, and thus the expectation 
values (Ek)- This is the sense in which quantum fluctu- 
ations are projected out, and it is this feature which we 
will restore in the following subsection. 

Comparing the kinetic energies for the uniform stripe 
configuration (Table IV, Fig. 11(a)), we find little varia- 
tion with J, and a strong preference for larger transverse 
hopping t y . Thus the effect of hopping anisotropy is as in 
the Hubbard model. Quantitatively, that 11% anisotropy 
leads to a 5% energy change, in contrast to the 30% effect 
in the Hubbard model with periodic BCs, can be ascribed 
to the fact that the Id state is already established in the 
isotropic case. The magnetic energies show a negligible 
dependence on the hopping, and a rather small overall ef- 
fect which is due to the large AF regions undisturbed by 
the presence of the stripe. However, the variation clearly 
indicates a preference for alignment of the stripe with the 
direction of stronger superexchange, a result readily ex- 
plained at zeroth order: the stripe cuts two bonds in the 
y-direction, and only one in x, so clearly costs less mag- 



netic energy in longitudinal alignment. Thus for filled 
stripes in the t-J model the two anisotropics compete. 
A detailed analysis by techniques-better able to account 
for quantum fluctuation effectstilO is clearly required to 
quantify this competition. 

Results for the half- filled, in-phase stripe (Fig. 11(b)) 
are given in Table V. The variation in kinetic energies 
with anisotropy is remarkably small in comparison with 
the uniform stripe, and it would appear fair to say that 
leading-order hopping contributions along the chain are 
cancelled by transverse hopping (primarily to the sites 
of larger hole density) for the ideal configuration estab- 
lished here. These results present no evidence to support 
alignment of the alternating stripe with the direction of 
strong hopping, and in fact show a weak preference for 
transverse alignment. The variation in magnetic ener- 
gies with anisotropy is similarly weak, as expected from 
the bond-counting argument. Corrections arise from the 
fact that the stripe now possesses a FM center line (c/. 
Fig. 9), and tend also to favor transverse alignment. The 
orientation of half-filled stripes is thus a delicate issue, 
and we are unable to make any firm statements on the 
basis of these results. For the same reason we have not 
tried to address further interesting questions which may 
be posed in this context, such as the energetics of corners 
in domain walls, crossings of stripes, or interactions be- 
tween stripes placed in proximity to each other .113 Finally, 
for both types of stripe we note that there is negligible 
interplay between the individual effects of anisotropics 
in t and J, with energy changes proceeding linearly and 
separately in both variables. 

B. Intermediate U 

In this subsection we will show results obtained in the 
spirit of setting U as large as possible while retaining the 
possibility of finding the true ground state. This is the 
approach adopted in Ref. |2^, and enables one to restore 
the competition between kinetic and magnetic energies in 
examining the above effects, but introduces a different set 
of questions concerning the physics of the t-U-J model.EZl 

We begin by characterizing the effects of J in the 
isotropic model with U = 50. Fig. 12 shows the evolution 
of the ground-state charge distribution as J is increased. 
Even rather small values of the superexchange interac- 
tion alter the J = structure (Fig. 12(a)) to one with 
stronger charge differentiation between and within small 
clusters (Fig. 12(b), also realized for J = 0.05). At in- 
termediate values of J there is evidence of competition 
with a phase (not shown) of small polarons showing pref- 
erential diagonal alignment. As J is increased towards 
the canonical cuprate-model value of 0.3, the small- J 
polaronic phase of Fig. 12(b) exhibits increasing spac- 
ing of the charge clusters, or preference for undisturbed 
AF regions (Fig. 12(c)), and shortly beyond this we 
find separated clumps of uniformly distributed charges 
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(c) (d) 

FIG. 12. Ground-state charge distribution for t-U- J model 
on a 12 x 12 cluster with open BCs, for hole doping x = 1/6, 
U = 50, and t x = t v = -1. (a) J x = J y = 0. (b) 
J x = J y = 0.1. (c) J x = J y = 0.3. (d) J x = J y = 0.5. 
Charge scale is such that largest circles represent 75.1% hole 
occupation. Spin distributions (not shown) are AF with un- 
altered phase. 

(Fig. 12(d)) in a phase which we regard as representa- 
tive of true phase separation. The quantitative value at 
which this separation occurs is not unexpected from pre- 
vious estimates 0~E3 

With periodic BCs this tendency to phase separation 
is strengthened, but so also is the tendency to forma- 
tion of anisotropic structures in the presence of hopping 
anisotropy (below). In certain parameter regimes it ap- 
pears possible for this anisotropy to suppress phase sepa- 
ration, which in a systematic study would lead to moving 
the phase boundary to higher J. By contrast, anisotropy 
in the superexchange has little effect at small J(< 0.1), 
where the system remains largely dominated by kinetic 
processes, but appears to promote the phase separation. 

In this parameter regime the spin state is not altered 
from its starting configuration of very weak antiferro- 
magnetism, which is simply reinforced during the self- 
consistency procedure. This maintenance of the spin con- 




(a) (b) 

FIG. 13. Ground-state charge distribution for t-U- J model 
on a 12 x 12 cluster with periodic BCs, for hole doping 
x = 1/6, U = 100, and J x = J y = 0.2. (a) tx = ty = -1- 
(b) t x = —1-1, t y = —0.9. Charge scale is such that largest 
circles represent 75.5% hole occupation. 

figuration is not an artefact of considering only a single 
spin component (Sec. Ill), as retaining all three com- 
ponents returns the identical result expected from spin 
rotation symmetry. While the resulting charge inhomo- 
geneities have in-phase nature, we cannot exclude the 
possibility that more extended structures (cf. Sec. IV), 
which reverse the spins over a finite region, form a sepa- 
rate class of minima inaccessible from this initial state. 

Fig. 13(a) illustrates this point with a result obtained 
in an isotropic system with periodic BCs and U = 100. 
The emergence of a diagonal stripe is a consequence of 
the commensurate filling, and in fact the structure should 
more correctly be described as a "bistripe" . The possi- 
bility offered by this filling to create a wall with a 2ir 
phase rotation (Fig. 14) allows this configuration to exist 
within an otherwise in-phase, AF background. Hopping 
anisotropy for this parameter choice leads to a rather 
different form of bistripe (Fig. 13(b)), whose orienta- 
tion is changed to vertical even by rather small differ- 
ences between t x and t y . As noted above, superexchange 
anisotropy J x ^ J y promotes polaron formation and 
phase separation of the above kind (Fig. 12(f)); quan- 
titatively, 11% anisotropy at J = 0.3 (cf. Fig. 12(c)) and 
22% at J = 0.1 (cf. Fig. 12(b)) are sufficient to drive the 
separation. Finally, the qualitative results described in 
this section depend rather little on the chosen value of 
the "large" U, which was varied between 16 and 100. 

We return now to the question of the competition be- 
tween the stripe alignment tendencies of anisotropics in 
t and J. The most unbiased approach we can find to 
investigate this is to set U = 16, and to begin with a 
uniform charge distribution, and zero initial spin on all 
sites except those in two corners. The following results 
are presented for a 12 x 11 cluster with periodic BCs, 
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FIG. 14. Ground-state spin configuration for t-U-J model 
on a 12 x 12 cluster with periodic BCs, for hole doping 
x = 1/6, (7 = 100, t x = t v = -1, and J x = J y = 0.2; 
corresponding charge distribution in Fig. 13(a). Spin scale as 
in Fig. 7. 



which we have chosen to favor stripes of the experimen- 
tal orientation. We used J = 0.5, as in Ref. [2^, but 
found little qualitative difference at J = 0.25 despite the 
results of Fig. 12. In Fig. 15(a) is shown the charge dis- 
tribution for parallel spins at diagonally opposite corners 
of the cluster, which for the chosen dimensions sets the 
requirement of an antiphase domain wall. Indeed a do- 
main wall is formed, but despite the commensurate filling 
(12 holes) it is not straight, and wandering of the domain 
line becomes more pronounced at smaller J, such that the 
line sections are closer to being diagonal. The competi- 
tion between a preference for diagonal domain walls from 
nearest-neighbor hopping and vertical domain walls from 
superexchange is exactly that expected from the simplest 
considerations on the idealized stripe. 

In Fig. 15(b) is shown the result of setting an in-phase 
configuration across the cluster: the holes prefer to form 
a diagonal domain line, and a wandering one whose local 
direction is again mostly diagonal. On periodic continu- 
ation of the cluster, the line appears as a sawtooth. The 
prevalence of kinetic terms over magnetic in the charge 
alignment is not apparent on the basis of Fig. 12, but 
we note that the doping here is much lower. In this case 
there is a reversal of the AF spin configuration across the 
domain wall, a result which argues for initial configura- 
tions containing a single "seed" spin, from which all other 
orientations may be "grown" . In fact such an approach 
produces few qualitative differences from the results pre- 
sented in this section, where the starting configuration 
was weakly AF. For smaller values of U the system is 
sufficiently flexible to choose the appropriate configura- 
tion, and for larger U the result of the "growth" proce- 
dure is either a rather more diffuse version of the same 
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FIG. 15. Ground-state charge distribution for t-U-J model 
on a 12 x 11 cluster with periodic BCs, for 12 holes (x — 1/11), 
U = 16, t x = t y = — 1, and J x — Jy — 0.5. Initial configura- 
tion with upper left and lower right spins (a) parallel and (b) 
opposite. Charge scale as in Fig. 3. Spin distributions (not 
shown) are AF with ir phase shifts across all domain walls. 



charge structure (with higher net energy), or a random 
arrangement indicating the local-minimum problem. 

Variations of the t- and J-anisotropies on the two 
states in Fig. 15 have little visible effect on the charge 
configurations even at the 22% level, although in 
Fig. 15(b) there is noticeable charge redistribution to a 
more "vertical" local stripe orientation. Thus we revert 
to analyzing the energies of the ground states, which are 
presented in Tables VI and VIII in the same format as 
in Tables IV and V. It is immediately clear that kinetic 
and magnetic energies are of the same order again in 
the regime around U — 16, and it becomes meaningful 
to consider their combination E to t- The next striking 
feature of the data in Table VI is the remarkably small 
change in the magnetic energy E m as a function of J x /J y : 
despite the fact that the system appears to contain a 
stripe, and has no frustrated bonds at the boundaries, 
variation of E m with J x / J y is in fact weaker than with 
t x /t y . Another interesting property is that the kinetic 
energy decreases for both types of hopping anisotropy, a 
result which would not be expected for a vertical stripe. 

The explanation for both of these features is contained 
in Table VII, which shows the components of Ek and E m 
for the x- and y-directions separately. The magnetic en- 
ergies (lower figures) vary linearly with J x and J y , with 
their sum (Table VI) almost a constant quantity. The 
kinetic energies change very strongly with t x and t y , as 
already observed in Table I, but again their sum changes 
only weakly. It is clear from the figures that the stripe in 
Fig. 15(a) is almost perfectly diagonal in character. The 
bond-counting argument for effects of J on stripe orien- 
tation give no preference in the case of a diagonal stripe, 
and show how well this picture works even for a system 
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TABLE VI. Kinetic (upper figure) and magnetic (lower) 
energy components per site for t-J model on a 12 x 11 cluster 
with 12 holes, periodic BCs, U = 16, and a starting configura- 
tion of parallel spins on opposite corners; charge distribution 
in Fig. 15(a). 





t x = -0.9 

ty = -1.1 


tx = -1 

ty = -1 


t x = -1.1 
t y = -0.9 


Jx = 0.45 
Jy = 0.55 


-0.56846 
-0.40620 


-0.55757 
-0.40098 


-0.56347 
-0.39812 


Jx = 0.5 

Jy = 0.5 


-0.56789 
-0.40712 


-0.55738 
-0.40160 


-0.56349 
-0.39848 


J x = 0.55 

Jy = 0.45 


-0.56733 
-0.40809 


-0.55718 
-0.40228 


-0.56350 
-0.39888 



TABLE VII. Energy components (E%), (E V K ), and 
(.Em) P er s ite for t-J model on a 12 x 11 cluster with 12 holes, 
periodic BCs, and U = 16. 





t x = -0.9 


tx = -1 


t x = -1.1 




ty = -1.1 


ty = -1 


ty = -0.9 




-0.20022 


-0.26859 


-0.33180 


J x = 0.45 


-0.36824 


-0.28898 


-0.23167 


Jy = 0.55 


-0.18584 


-0.18201 


-0.17987 




-0.22036 


-0.21896 


-0.21825 




-0.20045 


-0.26917 


-0.33278 


J x = 0.5 


-0.36745 


-0.28820 


-0.23070 


Jy = 0.5 


-0.20678 


-0.20255 


-0.20010 




-0.20034 


-0.19905 


-0.19838 




-0.20066 


-0.26975 


-0.33376 


J x = 0.55 


-0.36667 


-0.28743 


-0.22973 


Jy = 0.45 


-0.22778 


-0.22313 


-0.22037 




-0.18031 


-0.17914 


-0.17852 



with a broad charge distribution far from the idealized 
stripe. Changes in magnetic energy due to contributions 
from the hopping term are on the 1% level. The energy 
gain for both anisotropy directions also requires a pre- 
dominantly diagonal nature, although the minor prefer- 
ence for strong t y indicates that the stripe does preserve a 
small amount of (1,0) character (less than 5% compared 
to the changes in x- and y-components) . DMRG studies 
of the t-J model also show a significant presence of di- 
agonal stripes, and their near-degeneracy in energy with 
vertically oriented domain walls. till 

This result is rather surprising in view of the proper- 
ties observed earlier in the Hubbard model, and invite 
speculation on the relationship of t-U and t-U-J mod- 
els. A discussion of the energetics is incomplete without 
noting that on the same scale as in Tables VI- VIII the 
double-occupancy energy for U — 16 is approximately 
0.2 in every case, so is a quantity which may not be 
neglected. By contrast, in the previous subsection the 
double-occupancy energy for U = 1000 was on the order 
of 0.0036 in each case. 

Table VII shows the kinetic and magnetic energies un- 



TABLE VIII. Kinetic (upper figure) and magnetic (lower) 
energy components per site for t-J model on a 12 x 11 cluster 
with 12 holes, periodic BCs, U = 16, and a starting con- 
figuration of antiparallel spins on opposite corners; charge 
distribution in Fig. 15(b). 





t x = -0.9 

ty = -1.1 


tx = -1 

ty = -1 


tx = -1.1 

ty = -0.9 


J x = 0.45 
Jy = 0.55 


-0.55432 
-0.35716 


-0.56027 
-0.38007 


-0.56816 
-0.38864 


J x = 0.5 
Jy = 0.5 


-0.55126 
-0.35322 


-0.55928 
-0.37394 


-0.56619 
-0.38890 


J x = 0.55 

Jy = 0.45 


-0.55414 
-0.35327 


-0.55929 
-0.37393 


-0.56597 
-0.38927 



der anisotropy changes of 22% in t and J for the charge 
configuration of Fig. 15(b). Once again, the values are 
comparable, and they are also very close indeed to the 
figures in Table VI for a very different charge configura- 
tion; in fact it is the magnetic energy which would drive a 
preference for the broad, meandering stripe of Fig. 15(a). 
In Fig. 15(b) there appears to be a higher net alignment 
of the domain walls along the y-direction, and this may 
be taken to explain the preference shown by the kinetic 
energy for stronger t x . Here we notice a rather strong ef- 
fect of the hopping on the magnetic energy, for which the 
bond-counting argument also provides an explanation. In 
narrow stripes, a diagonal orientation breaks four bonds 
per site while a vertical one breaks three, and thus a 
hopping anisotropy (of either direction) which pushes a 
diagonal stripe towards a vertical one to gain transverse 
hopping energy will also lower the magnetic energy. That 
the net magnetic energy is lower for stronger J y is also 
consistent with the inferred (0,1) character. The fact 
that the numbers in the lower two rows match to very 
high accuracy is not an error, as shown by the fact that 
the separate x- and y-components of the relevant ener- 
gies (not presented here) have the same 20% and 40% 
differences as in Table VII, and rather reflects the deli- 
cate nature of the energy balance for the weakly oriented 
domain-wall structure. 

We have also considered the same calculations with 
different values of U{= 50) and of J{— 0.25), and find 
qualitatively similar results. Changing the filling x also 
leaves the general conformation unaltered, although the 
domain wall loops in Fig. 15(b) become shorter. The 
analysis for an 11 x 11 system returns the same spin con- 
figurations as in Ref. ^5], although we find the 12-hole 
domain line to be rather broad. For 12 x 12 arrays the 
domain lines tend to close, to avoid disturbing BCs which 
are already satisfied. 

We close this section with a summary of the robust 
conclusions which may be drawn from our studies of the 
t-J model. Without considering anisotropy, the superex- 
change term favors polarons over corrals. The presence 
of J also leads to a form of phase separation, namely the 
formation of charge clumps near the edges of an open 
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system. As a result, in the HF approximation, only 
systems with rather small values of J may show forma- 
tion of domain walls, whose alignment tends to be diag- 
onal. This tendency is particularly strong in the t-U-J 
model with intermediate values of U . Hopping anisotropy 
tx 7^ t y promotes uniform stripes normal to the strong 
hopping direction, as in the Hubbard model. There is a 
parameter regime of intermediate J which favors small 
polarons, but no clear sign of the non-uniform, in-phase 
stripes observed for the Hubbard model. Superexchange 
anisotropy J x ^ J y favors alignment of uniform stripes 
in the direction of stronger exchange, and thus competes 
with a hopping anisotropy of the same orientation. While 
quantitative results are hard to obtain, energetic consid- 
erations suggest that the relative effects of t anisotropy 
are stronger. For the boundary region of (t, J, x) parame- 
ter space between polaronic states and phase separation, 
we find that anisotropics in t and J act to stabilize co- 
herent structures. 



VI. DISCUSSION AND CONCLUSION 

In the previous sections we have detailed at length the 
many degrees of freedom allowed by the real-space HF 
approach, and the ways in which we use and interpret 
these to draw meaningful conclusions within this frame- 
work. We have also reviewed the experimental situation, 
focusing on those parts we are able to address, and have 
compared our results with existing theoretical literature. 
In this section we give a brief discussion of the more gen- 
eral context of our analysis. 

In Sec. II we noted that the LTT distortion of the 
LSCO system provides anisotropieS|-af 1-3% in the hop- 
ping and superexchange parameters. CJ Our investigations 
in Sees. IV and V required anisotropics on the order of 
10-20%, and these were not always sufficient to create a 
Id state. We appeal here to the quantitative insufficien- 
cies of the HF technique to argue that our results indeed 
suggest an important role for anisotropies, as suggested 
by the comparison of the hopping energy difference with 
other key energy scales. The superconducting T c , and the 
ordering temperature T c h of stripe formation, are both 
around 40K in LSCO. One of the very few indisputable 
results from the extensive studies of cuprate models is 
that there exist many competing candidate ground states 
whose energies are remarkably similar, and it seems not 
unreasonable that small anisotropies may play a signif- 
icant role in selecting among these. In addition to the 
more detailed energetic considerations and analytical or 
numerical approaches which we mention below, a more 
accurate assessment of anisotropy effects would also re- 
quire—moving from a one-band model to the three-band 
case.t£l This step is also a prerequisite for discussing at 
the microscopic level the issue of bond- or site-centered 
stripesfLrEZI which is one we have not mentioned hitherto. 
All of the charge structures we find by unrestricted HF 



for the one-band model, with and without anisotropy, are 
site-centered, and we do not regard the model as appro- 
priate for addressing this question. 

The connection between the lattice structure and 
anisotropies in the electronic system suggests the impor- 
tance of the phonon degrees of freedom. While dynamical 
phonons may play a significant role around the struc- 
tural phase transition, the static distortion requires a 
slow phonon mode which becomes soft. Although static, 
averaged phonon variables have been included in some 
HF analyses J23 we remain unaware of any studies where 
the coupled phonon modes which are introduced have 
allowed the possibility of a spontaneous anisotropy. As 
noted in Sec. II, we have studied only a model with 
fixed anisotropy set by the lattice, and considered the 
energetics of the electronic system. When the structural 
transition is in close proximity to the electronic one, the 
need for a coupled model becomes apparent. The rel- 
evance of lattice energies for the electronic system may 
also be clarified by detailed experimental analysis of the 
energetics of the structural transition. 

Other transition-metal oxide systems offer a clear ex- 
ample of the coupling of lattice and electronic degrees 
of freedom. The first observations of charge-ordering in 
stripe-like structures were made on nickelates, where the 
stripes are found to be diagonal (aligned along (1,1)). 
Detailed characterization of the ordering transition sug- 
gests once again a connection to the orthorhombic lat- 
tice structure arising from octahedral tilting at dopings 
x < 0.2.E3 Above this doping, charge ordering is a cou- 
pled structural and electronic event, which occurs at a 
single temperature, Tj, and this behavior persists to very 
high dopings (La2- K Sr x Ni04 with x < 0.5).E3ln this case 
the energy change of the electronic system on adopt- 
ing an anisotropic structure must be balanced directly 
against the distortion energy of the lattice, and theo- 
retical estimates of the competing contributions become 
more involved. In manganite systems, the most com- 
mon structural deformation is a Jahn- Teller distortion, 
which directly affects the electronic system by lifting the 
e g orbital degeneracyEd Such distortions usually alter- 
nate between sites ("AF"), and experimentally they ap- 
pear to influence the nature of the ground state primarily 
at low hole doping. However, observation of stripe-like 
structures in manganite systemsEiJ has not yet been re- 
lated unambiguously to lattice distortions, and may lie 
in the competition of next-neighbor repulsion V with 
anisotropic conduction-electron hopping t. 

With regard to methodology, the Hartree-Fock ap- 
proach is generally understood to be "crude but effec- 
tive" . It is capable of finding the relevant ground states of 
the system, and allows one to draw many qualitative con- 
clusions about their possible nature and order. However, 
the approximations involved in the breaking of symme- 
tries and neglect of quantum fluctuations are sufficiently 
strong that the parameter ranges where different features 
appear, i.e. the quantitative aspect, may not be trusted. 
In addition, with periodic BCs, the self-consistency pro- 
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FIG. 16. Schematic phase diagram of HF solutions for the 
Hubbard model as a function of U and hopping anisotropy et , 
for dopings 0.06 < x < 0.2. The "uniform" region is a para- 
magnetic metal. The nature of the stripe phases depends on 
x, as discussed in Sec. IV. Solid lines are true phase transi- 
tions, while the character of charge-inhomogeneous solutions 
at large U and small e t changes continuously. 

cedure has difficulties converging adequately due to the 
large choice of (translationally) degenerate ground states. 
A possible solution to both problems is, offered by the 
Configuration Interaction (CI) method.Ea L.3 In this ap- 
proach, the ground state is constructed as a (biased) 
linear combination of HF solutions. This technique, 
which has already been applied to certain aspects of the 
isotropic Hubbard model J^'Ell can be shown to give a 
systematic description of quantum fluctuation and tun- 
neling corrections to the HF solution, and of the quantum 
dynamics of the excitations. 

While CI represents a significant improvement to HF, 
there remain many other techniques which could be ap- 
plied to address the question of anisotropy effects at a 
more quantitative level. As listed in Sec. I, these include 
effective models of domain walls, the slave-boson repre- 
sentation, DMFT and DMRG, and we await with interest 
the results of subsequent analyses. We note finally that 
similar or adapted HF techniques may also be expanded 
to explore situations beyond the Hubbard and t-J mod- 
els. In the framework of more general t-J-LLV models 
oneanay seek finite, nearest-neighbor pairingcj and spin- 
flipO order parameters, with a view to discussing super- 
conductivity and flux phases in charge and spin. The 
latter possibilities have been of considerable recent inter- 
est in both single- and multi-band models. 

In conclusion, certain cuprate systems show a break- 
ing of 4-fold (x, y) symmetry in the Cu02 plane which 
is closely connected to the appearance of static stripes. 
Proceeding from this experimental motivation we have 
considered Hubbard and t-J models with hopping and 
superexchange anisotropy. The real-space Hartree-Fock 
approach offers many possibilities by which to gain in- 
sight into the nature of the solutions. In the isotropic 
Hubbard model we find predominantly closed, diagonal 



domain- wall loops (corrals), and polarons at larger U. 
Anisotropy quite generally brings a gain in electronic en- 
ergy, and introduces vertical stripes, which are stabilized 
by hopping between different charge densities. The qual- 
itative situation is summarized in Fig. 16, which without 
distinguishing between types of stripe is applicable for 
the doping range of physical interest. The dominant, 
uniformly "filled" stripes are oriented perpendicular to 
the direction of strong hopping, have antiphase nature, 
and show in their ground-state energy the very strong 
influence of transverse hopping. For certain hole densi- 
ties a non-uniform stripe appears, which may be the HF 
analogue of the "half- filled" stripes seen in experiment, 
and has longitudinal orientation and in-phase character. 
The isotropic t-J model in the HF approximation favors 
polaron formation and shows a tendency to phase separa- 
tion with increasing J. Anisotropics in t favor transverse 
stripe formation, while anisotropics in J reinforce longi- 
tudinal orientation of stable stripes, but also act to pro- 
mote phase separation. We conclude that stripes cannot 
be considered in isolation from the question of lattice- 
induced anisotropies. 
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